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(~| ■ Abstract. The definition of instantaneous eigenstate populations for a dynamical 

non-sclf-adjoint system is not obvious. The naive direct extension of the definition used 
for the self-adjoint case leads to inconsistencies; the resulting artifacts can induce a false 
inversion of population or a false adiabaticity. We show that the inconsistency can be 
^ I avoided by introducing geometric phases in another possible definition of populations. 

O^- An example is given which demonstrates both the anomalous effects and their removal 

by our approach. 
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1. Introduction 

^ I The adiabatic approximation is currently a commonly used tool to study quantum 

5^! dynamical systems p. Let s = t/T be the reduced time (T being the total duration 

which is assumed to be very large). Let \a,x{s)) be an instantaneous eigenvector of a 
selfadjoint hamiltonian H{x{s)) associated with a non degenerate eigenvalue Ea, i.e. 

H{x{s))\a,x{s)) = Ea\a,x{s)). (1) 

H depends on the reduced times s from some classical parameters x. Starting from the 
initial state iplO) = \a,x{0)), the state vector propagation obeys the time-dependent 



Schrodinger equation 

i dilj{sT) 



H{xismsT). (2) 



T ds 

The adiabatic theorem [H |2] states that the systems continuously follows the same 
eigenstate \a,x{s)) if the hamiltonian is gradually modified; that is 

^{sT) = Cais)\a, x{s)) + '^ci,{s)\b,x{s)) (3) 
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with 

Mh^a, - < 1, (4) 

|ca(s)p being called the population of the state a, if the following adiabatic criterion is 
fulfilled: 



Mh^a, |(6,a;(s)|^|a,a;(s))| 



{h,x{s)\ma,x{s)) 



« 1 (5) 



Ea{x{s)) - Ebixis)) 

which requires satisfaction of the gap condition: 3g > 0, Vs \Ea{x{s)) — Eb{x{s))\ > g, 
and also slow time variations within the hamiltonian H{x{s)). This leads to the 
wavef unction 

^{sT) ~ exp (^-ih-^T Ea{x{s'))ds'^ exp A{x{s'))ds'^ \a,x{s)), (6) 

where —ih^^T Ea{x{s))ds is the dynamic phase and — A{x{s))ds is the geometric 
phase. The geometric phase generator is defined as 

A{x) = {a,x\d\a,x) =^ A(x{s))ds = {a, x{s)\-^\a, x{s))ds. (7) 

Studies on adiabaticity and geometric phases can now be usefully extended 
to dissipative systems which use non- self adjoint hamiltonians to describe resonance 
phenomena [31 HI |5l [6l [7]. When resonances are defined by some scattering boundary 
conditions the resonance states are unbounded. This means that numerical and 
theoretical treatment of the resonances is a difficult task. Spectral deformation methods 
are commonly used to solve this problem. The first one, the complex scaling method 
[8], effectively transforms the resonance states into bound states. The hamiltonian 
becomes non-selfadjoint, and the resonances calculations produce non-real eigenvalues 
and the associated eigenvectors. The real part of a resonance eigenvalue corresponds to 
the resonance energy and its imaginary part corresponds to the resonance width (the 
inverse of the resonance lifetime). A second approach is the optical potential method [9], 
which has similar properties. These two methods of resonance modelling involve non- 
selfadjoint hamiltonians and are used successfully to treat molecular photodissociation 
problems [TOl |TTl [12]. The present paper deals with non-hermitian matrices. Such 
matrices can be viewed as effective hamiltonians associated with the hamiltonians 
of resonance problems. Effective hamiltonian techniques use small hamiltonians 
to describe complex systems without loss of informations. The principal effective 
hamiltonians techniques are the partitioning technique [13], the quantum KAM method 
|14j . the adiabatic elimination method [H] and the Bloch wave operator method [13] . 
Small non-selfadjoint matrices have also been used as hamiltonians to describe some 
photoionization phenomena modelled by bound states coupled with a structureless 
continuum (these states are then viewed as resonances) [151 [HI EZ] ■ Finally we should 
note that the formalism used in this paper can be used to treat open quantum systems 
[TS] [T9] which are described by a Lindblad equation 
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where p is a density matrix (a traceless positive selfadjoint matrix) and the Lindbladian 
L is a "superoperator" . For a finite dimensional Hilbert space, dim'H = n, the 
space of the density matrices can be identified with an n^-dimensional Hilbert space 
£ (usually called the Liouville space). In the Liouville space, the Lindblad equation 
takes the form of the usual Schrodinger equation. If the Lindbladian is reduced to 
be a commutator with a selfadjoint hamiltonian, i.e. L{p) = [H,p\ then L is a 
selfadjoint operator in the Liouville space C. Usually, however, the Lindbladian is 
L{p) = [H,p] — ^J2ki{^k^k,p} + SFfcprj,) where the operators model the different 
decohering processes ({., .} is the anticommutator). In this case, L is a non- selfadjoint 
operator in C. The non-real eigenvalues of L (its resonances) are associated with the 
decoherence. The present work can be applied to a such system, but for the sake of 
simplicity we will use only the name "hamiltonian" even though it could represent a 
Lindbladian in the Liouville space. 

For non-self adjoint hamiltonians, the adiabatic theorem [201 HI] involves a criterion 
on the dissipation rate of the quantum system studied (i.e. on the imaginary 
part of the instantaneous eigenvalue considered). Naively, we could imagine that 
the adiabatic approximation for dissipative quantum systems should take the form 
ipisT) = J2b'^b{s)\b, x{s)) with \fb ^ a, ^ ^1. The following section shows that 
this approximation is inconsistent, since the populations |cbp are not well defined. 
Consequently a better definition of the instantaneous populations is proposed in 
section [21 This definition is equivalent to the c-product condition (described in [3]) 
for the case of symmetric non-hermitian hamiltonians and remains valid even for the 
non-symmetric case. The illustrative example of a 2-state dissipative system is then 
described in section [31 with an application to an adiabatic state fiip [221 [23] which is 
generated by following a loop around an exceptional point in the parameter plane. 



2. Populations of the instantaneous eigenvectors 

For the sake of simplicity we consider a non-selfadjoint hamiltonian H{x) of rank 2 (the 
discussion can easily be generalized to higher dimensions). We deal with a biorthogonal 
basis set and we use the standard scalar product. Let {Ei{x) , E2{x)} be the two 
instantaneous eigenvalues of H{x) (assumed to be diagonalizable) and {|1, x), |2, x)} 
the two associated eigenvectors, and let \2*,x)} be the biorthogonal basis, that 

is: 

H{x)\a,x) = Ea{x)\a,x), (8) 
H{xy\a*,x) = Ea{x)\a*,x), (9) 

{a*,x\b,x) = 6ab- (10) 

(The bar denotes the complex conjugate). We note that in certain highly symmetric 
cases (i.e. when if* = H) the {|a, x)} and the {|a*,x)} are related by a simple complex 
conjugation rule (formalized by the c-product [2^ [3]). 
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Let 

^ ' \ (2*,x|d|l,x) (2*,x|d|2,x) ' ^ ' 

be the matrix of the geometric phase generators and of the non-adiabatic couphng. Let 

VXt) = ?7(t,0)|l,x(s = 0)) (12) 

be the wave function for an evolution s i— )■ x{s) {U{t,0) being the evolution operator). 
As in conservative systems, we could write 

V^(sT) = ci{s)\l,x{s)) + C2{s)\2,x{s)), (13) 

|ci(s)p being the "population" of the state 1 and |c2(s)p the "population" of 
the state 2; the adiabatic approximation corresponds to the case where |c2p is 
negligible. However this definition presents a significant problem. By constrast with 
the conservative case, where we choose an orthonormal eigenbasis to represent the 
dynamics, the normalization convention is arbitrary for a biorthogonal basis. Indeed, if 
({|l,a;), |2,a;)}, {|l*,a;), |2*,a;)}) constitutes a biorthogonal eigenvectors system, this is 
also the case for (^{\l,x),\2,x)},{\l*,x),\2*,x)}^ with 

\a,x) ^Xa(x)\a,x) (14) 
\a*,x) ^ {Xa{x))~'^\a*,x) (15) 
for all arbitrary functions Xaix) G C*. We must then set 

ij{sT) = gi(s)|M5)) + C2{s)\2^)) (16) 

with 

The population |cap = j^j- then depends on the arbitrary normalization convention. 
The calculation of the instantaneous populations then makes no sense, since the 
adiabatic multiplier ci can artificially grow if the normalization of \a, s) decreases. 

To solve this problem, wc introduce a definition of the instantaneous populations 
using geometric phases. We set 

4;{sT) = di{s) exp ^ Au{xis'))ds'^ |1, x{s)) 

+ d2{s) exp ^- ^ A22{x{s'))ds'^ |2, x{s)) (18) 

with {a,x{s = 0)\a,x{s = 0)) = 1. The instantaneous populations are now defined by 
|da(s)p. Immediately we notice that for the conservative case this definition coincides 
with the standard definition of the instantaneous populations since |e~ -^o" ^aa(a;(s'))ds' |2 _ ^ 
in a self-adjoint system {Aaa is purely imaginary; the standard definition of the 
instantaneous populations is invariant under arbitrary changes to the phase convention 
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of the eigenvectors). By making an arbitrary change in the normahzation convention, 
we have 

Aaa{x) ^ {0'*,x\^\a,x) (19) 

= , , / + a*, a; d a, x) 20 

= dlnA„(x) + ^„„(x). (21) 

We then have 

g- Aaa{x{s'))As' ^ K{X{())) ^_ js A„„(a;(s'))d.' _ /22) 

Aa(a;(s)) 

In order to preserve the initial condition, we set Xa{x{0)) — 1. This leads to 

iP{sT) = di(s)e--''o ^iiW-'))d.'|{^)^ + d2(s)e-^o ^22(x(s'))d.'|2'^)^_ (23) 

The definition of the instantaneous populations, |(iap, is now invariant even in the event 
of arbitrary changes in the normalization convention. Here we insist on the fact that 
the coefficients da are more intrinsic than the coefficients c^. 

Owing to dissipation the total population strays away from 1 (except for s — 0). We 
have 1 — (|(ii|^ + |(i2n > 0, but this is not a properly defined dissipation rate (i.e. 
\di\^ + |o?2p 7^ llV'in- The dissipation rate is given by (1 — ||'0(sT)|p) with 

msT)f= V4R4(s)expf- [\Aaaix{s')) + A,,{x{s')))ds' 

a,b ^ 

X {b,x{s)\a,x{s)) (24) 
= ^db{s)rjba{s)da{s) (25) 



a, 6 



D{sy7]{s)D{s) (26) 



di(s) 



with D{s) = ^ and r]ba{s) = c-/o(^-(^(*'))+^«(^(^')))d^'(6,a;(s)|a,a;(s)). rj 

\ d'2[s) J 

constitutes an s-depcndcnt scalar product for D{s), the representation of the wave 
functions in the instantaneous (non-orthonormal) eigenbasis. We note that rj is well 
defined since it is independent of the normalization convention, fj = rj. The scalar 
product matrix satisfies the following differential equation: 

^^A^rj + rjA (27) 
ds 

with 

Aab{s) = (a*,a;(s)|e^o^-W«'))d.'^ (^e-io'^'"'(^(^'»^^>,a;(s))) . (28) 

We can then write 

r]{s) = (Te^oMs')dsy ^^Q^j^f^A(s')ds' (29) 

where Te is the s-ordered exponential (the Dyson series). 

Finally, if we want the instantaneous total population to be consistent with the 
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dissipation rate, then the instantaneous population must be defined as a{s)\da{s)\'^ 
with a = -^Y^. 



Using this analysis, we propose that the consistent adiabatic approximation is 
d2 



di 



<C 1 



(30) 



and not 



<^ 1. In a similar manner, the adiabatic criterion must be independent of 



the arbitrary normalization convention: 



^J^ 5Je(A,,(.(s'))-A.„(.(s')))d.' I ^^^^ I ^ |a, xis))\ 



mAt,t{x{s'))-Aaa{x{s')))ds' 



< 1. 



{b*,x{s)\^\a,x{s)) 



Eaix{s)) - E,{x{s)) 



(31) 



We should point out that the use of the parallel transport condition, i.e. the 
definition of the populations with eigenvectors {\a,x{s))}a and {\a*,x{s))}a such that 



{a*,x{s)\ — \a,x{s)) = 



(32) 



implicitely takes into account the instantaneous population definition which includes 
the geometric phases. Indeed, it is easy to show that the {\a,x{s))}a are related to an 
arbitrary set of eigenvectors {\a,x)}a by 

= e-^o^'^^(-i^'^)'^^'\a,x{s)). (33) 

In the particular case of a symmetric matrix, the left eigenvector are the complex 
conjugates of the right eigenvectors, so that the parallel transport condition is also 
equivalent to the c-product normalization condition for the eigenvectors |a, x(s)c.p.), 
fixed by [3] 

(a,x(s)c.p.|a,x(s)c.p.) = 1. (34) 

The result of this paper could be expressed without explicit reference to the geometric 
phases by saying that the treatment of the population tracking for non-hermitian 
systems necessary needs to impose the parallel transport condition (whereas this is not 
necessary for hermitian systems). Nevertheless, we prefer here to make the geometric 
phases appear explicitely. Indeed the {\a, x{s))}a are defined only along the path in the 
parameter space defined by s i— )■ x{s)] for a different path s i— )■ x'{s) the eigenvectors 
{|a, x{s))'}a are different. It is then more general to consider {|a, x)}a (without a parallel 
transport condition) which are defined globally on the whole of the parameter space. 
Moreover we note that we cannot use the parallel transport condition to define the 
eigenvectors if the path is closed x{s = 0) = x{s = 1), because of the double definition 
of the eigenvectors at a;(0) = a;(l): 

lo^^)) =e-^o^-(-W)d«|a,x(l)) ^ |a,x(0)) = |a7^)). (35) 
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The impossibility of giving a single definition of the eigenvectors in the parallel trans- 
port condition e~ Aaa{x{s))As jg Qg^Hg^j holonomy of the parallel transport. 

The next section illustrates that numerical artifacts due to bad definitions of the 
populations can induce a false adiabaticity or a false inversion of population in numerical 
simulations. 

3. Illustrative example 

The definition (fT8l) is relevant for non-hermitian symmetric or non-symmetric matrices. 
While the formalism used here should work for larger more realistic hamiltonians and 
for Floquet-type hamiltonians, we study here a low dimensional illustrative example for 
which the data can be shown conveniently. Let the parameter dependent hamiltonian 
be 

{w,z) G C^. This hamiltonian is associated with a quantum bound state coupled to 
a quantum resonance (with resonance width F) by a laser field with amplitude VL and 
phase 0. The laser field is quasiresonant with the transition from the bound state to 
the resonance with a detuning value equal to A. We assume that we can modulate the 
complex numbers {w, z) to generate a dynamics. 



3.1. False adiabaticity 
The spectrum of H is 

Ei{w, z) = z — + z"^ = z — V (37) 

E2{w, z) = z + v^l^P + = z + v (38) 



where v = ^y\w\'^ + z'^ = zyl + (we choose the Riemann sheet such that = z). 
We restrict our attention to the parameters {w, z) such that Qm{E2 — Ei) = 2Qm{v) < 
(£"1 is the less dissipative) . 

We can easily verify that the eigenvectors of H{w, z) are 

\l,w,z)=^^(^'^J^ (39) 

|2,t.,^)=72(^^_~J j (40) 
where 7a are the appropriate factors to fix the initial norm of the basis vectors to 1, i.e. 

71 = ^ (41) 

V(k(0)+t;(0)|2 + !«;(— 
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72 = ^ (42) 

vTRoT^WTRo)?) 

The biorthogonal basis set is 

1 

2v{v + z)^l 



l*,w,z) = ^ (43) 



\2*,w,z)^—^^-^^\2,w,z). (44) 
2v(v - z)7j 

The associated geometric phase generators are 

iBw + WW WW {z-\-v)z , , 

^ Av^ + 2v{v + z) + ^^^^ 
WW + WW WW {v — z)z 

^''^ Av^ ^ 2v{v-z)~^h^- ^^^^ 



We can verify that another possible eigenstate 1 reads 

7i / w 
/3 [ z-v 



17^) = ^ ^ I (47) 



with 



and 



\l,w,z) = P^^^ll^z) (48) 



w 



w{0) 



v{0) + z{0) 



(49) 



\l*,w,z) ^ ^^^J^' 2 \'^,w,z). (50) 



The factor P ensures that 11,^(0), z{0)) — \l,w{0), z{0)) in order to preserve the initial 
condition. We assume that /5 7^ 0. The associated biorthogonal eigenvector is 

2v{v — z)^l 

The normalization factor between the two conventions approaches zero at the limit 

\w\ 0, 

w 

lim (5 = 0. (51) 

\w\^0 V + Z 

With this normalization, the first geometric phase generator becomes 

J WW + WW WW {v — z)z . . 

^''^ ^ 2v{v-z)~^^- ^^^^ 

These properties induce the following analysis. Let 

iP{sT) = di{s)e-^o^'^(^'}<^''\l,n^z{s)) + d2is)e-^o^''^^')^^^ (53) 

be the solution to the Schrodinger equation for the evolution s i-> {w{s), z{s)) with the 
initial state ijj{0) — \l,w{0), z{0)) . We suppose that the adiabatic approximation is 
^1. If we consider the naive definition of the population by setting 

^|;{sT) = c,{s) \ 1, w{s), z{s)) + 02(5) |2, w{s), z{s)) 



valid, i.e. 
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then the quotient 



Cl 


1 \w\ 


Cl 


C2 


(5\z + v\ 


C2 




1 \w\ 


di 




P\z + v\ 


d2 



Jo^eiA22is')~Aii(s'))ds' 



(54) 

(55) 



with ip{sT) = ci{s)\l,w 
5Re(A22 - Au) = and 

Iwl < 1 = 



S),z{s))+C2 

— I'"! 



di 



m\z+v\ 

^ 1 and 



2,w{s),z{s)). We note that if Q^m(w) = then 
. The hmit f lSTj) can induce 



< 1. 



(56) 



Thus we observe a false non-adiabaticity due to the badly defined populations. In this 
case we also have a false population inversion (in the sense that population 1 is negligible 
with respect to population 2 using the badly defined populations, while the well-defined 
populations give the inverse result). 



Conversely, let ip{sT) [cf. ( 153|) ] be the solution to the Schrodinger equation for an 
evolution starting with an initial state '?/'(0) = |2, t(;(0), 2;(0)) such that the adiabatic 

for s sufficiently large, \di\ ~ \d2\. By the same 



approximation is not satisfied, i.e. 
arguments, if |w| ^ 1 we have 



di 
d2 



1 and 



< 1. 



(57) 



In this case, we obtain a false adiabaticity due to the badly defined populations. 




Figure 1. Time variations of the Rabi frequency w (line) and of the detuning 5fte(z) 
(dashes) . 



Figure [T] shows a simple example of a gaussian variation for w{s) (coupling) and a 
decreasing cosine function for 9?e(z) (detuning), so that \w{t)\ ~ near the end of the 
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Figure 2. Top: |ci(s)| (line), |c2(s)| (long dashes), |di(s)| (short dashes), |d2(s)| 
(dots); bottom: ratios |c2(s)/ci(s)| (line) and \d2{s) / di{s)\ (dashes). The initial state 
is |l,u;(0),z(0)). r = 0.1. 



1 

0.8 
0.6 
0.4 
0.2 


0.2 0.4 0.6 0.8 1 

1 

0.8 
0.6 
0.4 
0.2 


0.2 0.4 0.6 0.8 1 

Reduced time s 

Figure 3. Top: (1, w, z|l, u>, z) (line), (1, w, z|l, u>, z) x 

|g/o=Aii(s')ds'|2 (dashed line). 

Bottom: (2, w, z|2, w, z) (line), (2, w, z|2, z) x |e/o ^22(s')ds' |2 (dashed line). This 
corresponds to the case of a fake inversion of population with F = 0.1. 





considered time interval [0, 100]. The chosen functions are 

w{s) = Woe 2^ (58) 
A(s) = ^e{z{s)) = Ao cos(0.47rs) (59) 

with wq = 1, Aq = 0.5, r = 0.1 or 0.2 and a = 0.16. This elementary example 
perfectly illustrates our assertions. Figure [2] corresponds to the case of a false inversion 
of population, beginning with the state \l,w(0), z(0)). 1 — 1 becomes very large when 

^1 

s > 0.8, although |^| remains very small for the duration of the interaction (using 
the definition which takes the geometric phases into account to compensate for the 
unstable norm of the basis vectors). This is consistent with our analysis because this 
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Figure 4. Top: |ci(s)| (line), |c2(s)| (long dashes), |di(s)| (short dashes), |d2(s)| 
(dots); bottom: ratios |ci(s)/c2(s)| (line) and \di{s) / d2{s)\ (dashes). The initial state 
was |2,u;(0),z(0)). T = 0.2. 



false inversion of populations occurs when w{s) becomes very small due to the gaussian 
function. Moreover we note the uncontrolled increase in ci and C2, both reaching values 
larger than 1. 

Figure [3] shows the principal cause of this problem (the uncontrollable variations 
in the norms of \l,w,z) and \2,w,z)) and also the exact compensation obtained with 
the exponential of the geometric phases, leading to two stable unitary norms (since we 
have taken the precaution of setting the initial norms to one). 

In contrast, figure H] illustrates the inverse phenomenon of false adiabaticity. The 
initial state is \2,w{0), z{0)). Thus the ratio | — | stays under the value of 0.01 as if it 
were an adiabatic evolution, while | ^ | increases, with two components of the same order 
of magnitude {di ~ ^2) at the end of the interval. 

3.2. Closed loop around an exceptional point in the parameter space 

If r and are fixed, the matrix of ( l36l) has two exceptional points in the parameter plane 
(n. A). They are located at (fi, A) = (±r/4, 0). At these points the two eigenvalues 
and the two eigenvectors coalesce. If a closed loop that encircles a single exceptional 
point is followed in the parameter space, then the adiabatic basis is transported so as 
to have the following property [1]: 

|l(s = l)) = i^i|2(s = 0)) 

|2(s = l)) = z/2|l(s = 0)) (60) 

where ui and z/2 are complex numbers. If the dynamics is adiabatic, starting with one 
of the eigenstates |l(s = 0)) leads to a final state |l(s = 1)) which is proportional 
to |2(s = 0)). This interchange is called an adiabatic flip. This interesting form of 
degeneracy in non-hermitian quantum problems is the subject of several recent papers. 
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involving both theoretical [221 ESI ESI ESI EZ] and experimental ^ EH] studies. The 
references [221 ES] deal with a symmetric matrix and show that beginning the loop 
with one state is favourable to an adiabatic behaviour which leads to an adiabatic flip 
(interchange), while beginning with the other one (the more dissipative one) induces 
strong non-adiabatic couplings and state exchange during the process, whatever the 
time duration. The occupancy coefficients of the two states and the states themselves 
are exchanged during the loop, leading finally to an absence of fiip. This is consistent 
with the conditions associated with the applicability of the adiabatic theorem [201 El]- 
We can now confirm these results by our calculations and also extend them to the 
non-symmetric case by using the definition (ITSll . 

3.2.1. Symmetric case We choose the parameters 

T = 100, 

z{s) = A - ir/4 

with r = 0.5 and A(s) = 0.24 T sin(27rs), 

w(s) =fi(s) = ^ + 0.24rcos(27rs). (61) 

We calculate the different adiabatic multipliers Ci(s), C2(s) and di{s), d2{s) following the 
adiabatic eigenstates by using the two conventions defined above. We also calculated 
results using a third convention which corresponds to the c-product normalization (IM]) 
and which gives the coefficients 

ei = {l,wis),z{s),.^,\i;{sT)), (62) 

62 = {l,w{s),z{s),.p,\tlj{sT)) (63) 

such that 

\^{sT)) = e,{s)\l,w{s),z{s),,p) + e2{s)\2,w{s),z{s),,p,). (64) 

This c-product normalization is adjusted so as to begin with the same initial condition 
(ci, 62) = (1, 0) (or (ei, 62) = (0, 1)) as the two others. Figures [5] and [6] correspond to the 
dynamics issuing from the initial states \l,w{0), z{0)) and |2, w(0), ^(0)), respectively. 
The |ci| and \c2\ curves seem to indicate a non-adiabatic behaviour but are not 
numerically significant. The evolution of ei and 62 is exactly the same as the evolution 
of di and d2- We see an adiabatic evolution in figure [5] for the dynamics issuing from 
\l,w{0), z{0)) . The final state |1, w(l), ^(1)) is predominantly occupied at the end of 
the loop but owing to the fiip of the eigenstates we do have an adiabatic fiip. When the 
initial state is |2, w{0), z{0)) (figure[6l), the non-adiabatic exchange of populations during 
the first part of the loop induces a final state near to state |1, w{l), z{l)) oc |2, w{0), z{0)) 
and there is no fiip. 

3.2.2. N on- symmetric case The adiabatic loop is the same but the only difference from 
the previous case is that the off-diagonal elements are now complex conjugates, so that 
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Figure 5. Top: (line), |(i2(s)| (dashes); middle: |ei(s)| (line), |e2(s)| (dashes); 

bottom: |ci(s)| (line), |c2(s)| (dashes). The initial state is |1, w(0), z(0)). 
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Figure 6. Same as figure[5]with initial state |2, w(0), z(0)). 



the left eigenvectors are not the complex conjugate of the right eigenvectors. We set 

w{s) = Q{s) e'^ = + 0-24 T cos(27rs)^ e'^ (65) 

Figures [7] and [S] show the evolution of the occupancy coefficients with the different 
definitions, beginning respectively with the initial state |1, w{0), z{0)) and |2, w{0), z{0)). 
We clearly see that it is no more possible to use a c-product-type normalization: fl^ 
should not be used in this case. The coefficients Ci and C2 continue not to be numerically 
significant. Only the coefficients di and d2 can be used in the present case and they 
show the same behaviour as for the symmetric case. 
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Figure 7. Top: (line), |(i2(s)| (dashes); middle: |ei(s)| (line), |e2(s)| (dashes); 

bottom: |ci(s)| (line), |c2(s)| (dashes). The initial state is |1, w(0), z(0)). 
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Figure 8. Same as figure[7]with initial state |2, w(0), z(0)). 



4. Conclusion 

This paper should be regarded as a simple clarification about the difficult interpretation 
of calculations on dissipative quantum systems described by non-hermitian hamiltonians 
and about the use of state populations as important observables, especially for studies 
on adiabatic phenomena. When one works with dissipative quantum systems, the 
erratic time evolution of the norm of the adiabatic eigenvectors can create difficulties 
in answering the question: what is the relevant definition of a population? When 
the eigenvectors are calculated numerically (for large matrices), there is no reason for 
them to be continuously transformed from one point to another in time, thus the norm 
variations can appear to be quite disordered (worse than in our semi-analytic example). 
Seemingly, the easiest intuitive solution is to artificialy normalize the "right" eigen- 
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vectors as if we were working with an orthogonal basis set, but this is not a coherent 
way to work. We have shown that it is much preferable to compensate for the erratic 
variations in the norm of the adiabatic basis set by including the exponentials of the 
geometric phases in the basis vector decomposition, leading to an invariant definition of 
the populations under arbitrary changes of the normalization choice. 
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